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Abstract. In this paper, we present several gravity tests made in using the last 
INPOP08 planetary ephemerides. We first propose two methods to estimate the 
PPN parameter /3 and its correlated value, the Sun J2 and we discuss the correla- 
tion between the Sun J2 and the mass of the asteroid ring. We estimate possible 
advance in the planet perihelia. In the end we show that no constant acceleration 
larger than 1/4 the Pioneer anomaly can affect the planets of our solar system. 

1. Introduction 

Thanks to the high precision achieved with the observations deduced from the tracking of 
spacecrafts, it becomes possible to estimate relativistic parameters, for instance mainly 
7 and /?, of the Parametrized Post-Newtonian formalism of General Relativity (Will, 
1993). Nevertheless if 7 plays a role in the equations of motion, it is worth to note that 
light propagation is only sensible to that parameter. PPN 7 can be then estimated with 
high accuracy by light deflection measurements by VLBI (Shapiro et al. 2004), by time 
delay during an interplanetary roundtrip and by Doppler tracking data of a space mission 
(see for instance the Cassini experiment, Bertotti et al. 2003). This is also why, in the 
following, we put 7 = 1 in order to only test the sensitivity of PPN (3 on the perihelion's 
advance of planets. However a physical quantity relative to the Sun, its oblateness J2, 
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plays also a key role in this phenomena. Indeed, usual expression of the advance of 
perihelion is given by (Will 2006) 

a(l-e2)c2 a2(i_g2)2 ^ > 

where G and c are the newtonian gravitational constant and the speed of light in a vac- 
uum, respectively; J2 and Rsun ai'c the Sun oblateness and the Sun equatorial radius; 
a and e are the semi- major axis and the eccentricity of the precessing planet. One can 
immediately suspect that it is not possible to do a relevant estimation of PPN (3 with- 
out considering Sun J2. Furthermore, it may be not possible to decorrelate safely these 
two quantities with only one planet. With INPOP08 (Fienga et al. 2009), MEX and 
VEX tracking data have lead to an important improvement of Mars and Venus orbits, 
respectively. It is then suitable to take advantage of this new situation by attempting to 
decorrelate these parameters. 

Indeed, the impact of the Mars and VEX observations is not only limited to the 
improvement of the planet dynamics. They also play a role in the determination of pa- 
rameters such as the asteroid masses, the oblateness of the Sun and the PPN parameter 
(3. The ratio between the uncertainties of the observations and the sensitivity of the 
observed orbit to the GR modifications has been evaluated by dividing the cumulative 
advance of the perihelion over a period of time corresponding to the time span of observa- 
tions by the angle uncertainty of INPOP and presented in table[TJ If the amplitude of the 
advance of the perihelion on Venus and Mars orbits is considered for a set of observations 
of equivalent accuracy, Venus data will be seven times more efficient to test general rela- 
tivity and to estimate the sun J2 than Mars. If VEX mission is prolongated from 2 years 
to 4 years and if VLBI observations are done from the tracking of the spacecraft with an 
accuracy of about 1 mas, VEX data will be then as important for the PPN testing and 
Sun J2 estimations as the direct 800-meter accuracy radar ranging on Mercury. Besides, 
the Mars data are still very important because of the long time span of observations of 
very good quality obtained since the Viking mission in 1978. 

Thanks to the informations brought by the combination of very accurate tracking 
data of spacecraft orbiting different planets, the planetary ephemerides become then an 
interesting tool for gravity testing. In the following, we give some examples of such tests. 



2. Determination of PPN /? and sun oblateness J2 

2.1. Correlation between sun J2 and asteroid modeling 

The advance of the perihelion induced by general relativity and sun J 2 has an impact 
very similar to the advance induced by the main-belt asteroids on inner planet orbit. In 
INPOP08, a ring was fixed to average the perturbations induced by the main-belt aster- 
oids which cannot have their signal fitted individually on tracking observations. This ring 
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Table 1. The first 2 columns give the a-priori INPOP uncertainties in geocentric angles 
and distances limited by the observation accuracies. In the third column, one may find the 
estimation of the general relativity and sun oblateness effect of the advance in perihelion, 
uj, on the Mercury, Venus and Mars perihelia per year. The fourth column gives the S/N 
ratio estimated over the period of time given in years in column 5. 
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has its physical characteristics (mass and distance to the sun) estimated independantely 
from the fit by considering the albedos and physical properties of 24635 asteroids (for 
more details see for instance Kuchynka et al. 2008). 

As illustrated on figure [l] there is a correlation between the effect on the geocentric 
distance of the modeling of the ring as done in INPOP08 in one hand and the effect of the 
sun oblatness in the other hand. Indeed, on these plots, one may see how a small change 
in the value of the sun J2 (12%) induces after the refit of the planet initial conditions 
a periodic effect very similar in amplitude and frequency on Mercury, Mars and Venus 
distances to the Earth than a change in the mass of the asteroid ring (17%). Besides, 
the Saturn-Earth distances are not affected in the same way. As the ring does not affect 
the outer planets as it does on the inner planet, it becomes possible to decorrelate the 
signal induced by a small change in J2 from the one procuced by the ring when one can 
note an inversion of the GM signal for example in 2006. This result is consistent with the 
analytical study done by lorio (2007) concluding also to the importance of taking into 
account the effects of asteroids on planetary orbit during relativistic tests. 

It stress also the crucial importance of having a modeling of the asteroid perturbations 
as a fixed ring characterized independantly from the fit of planetary ephemerides. 

We limit then a surestimation of the value of the sun J2 melting in this value some 
effects induced by the asteroids. 

Two different but complementary analysis and determination of PPN f3 and Sun J2 
are presented in the next sections with a fixed model of asteroid perturbations (same 
values of asteroid and ring masses and of densities as INPOP08). 

2.2. Estimations by least squares 

The first approach is based on the classic least square estimation of parameters during 
the fit of planet equations of motion to observations. To check numerically the simplied 
assertion made in introdcution (section [1]), we estimate here what is the impact of each 



4 



Fienga et al: Gravity tests with INPOP 



• Delta_GM_ring 



8 10 



Delta J2 
Delta_GM_rjng 



Earth-Venus 



8 10 



Delta, J2 
Delta_GM_ring 




Delta J2 
Delta_GM_rjng 



8 10 



Fig. 1. Residuals obtained by comparisons between MGS/MO, VEX and Cassini range 
tracking data and ephemerides perturbed by a small change in sun J2 (12%) and by a 
small change in the mass of the asteroid ring (17%). 



datasets in the determination of J2 and (3: several adjustements of the initial condition 
of planets and the parameters J2 and /3 are made using different sets of observations. 
This leads to 32 adjustments based on INPOP08. For each fit, changes were made in the 
selection of Mars and Venus data in order to estimate the impact of each important set 
of observations in the fit of the Sun J2 and PPN (3. We look at the variations in the 
estimation errors of the 2 parameters and we use the 1-cr given by the least squares as 
indicator of this uncertainty. With this method, we are then able to quantify the influence 
of each data sets on the determination of the pair (/3, J2) as well as the stability of the 
determinations of the parameters. Indeed these variations in the error's estimation of the 
pair (/3, J2) are a relevant indicator of the uncertainty of the fit of [3 and J2. 

To take into account the correlation between J2 and [3, we use two modes of adjust- 
ments: in the model, (3 or J2 are fitted alone with the initial conditions of planets; in the 
mode2, both parameters are fitted simultaneously with the initial conditions of planets. 

The results are summarized in table [51 One can first notice that the determinations of 
sun J2 and /? made separately (i. e. model) give better a than fits including simultane- 
ous (/3, J2) determination (mode2). This is obviously consistent with the expected result 
relative to the determination of correlated parameters. The best results for a correlated 
determination of J2 and /3 (mode2) are then obtained when only the most accurate obser- 
vations of Mars (MGS/MO, MEX and Vking) and Venus (VEX) are used simultaneously. 

Moreover, we note that the combined use of Venus ranging data and the complete 
data set for Mars do not really improve the separated determination (model) of /? and 
J2, mainly due to the low accuracy of these observations, but a contrario it gives better 
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Table 2. 1-a least squares obtained for 32 and (3 using several sets of observations. 
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correlated estimations (mode2). This is also consistent with the fact that fitting over 
observations from two different planets helps to decorrelate safely J2 and f3. Furthermore, 
the Viking data by prolonging the fit interval with observations of rather good accuracies 
allow a decreasing of the uncertainties of about 20 % for J2 and about 10% for (3. 
Finally, it appears that the VEX data improves the determinations in a significant way: 
decreasing of 31 % of the least squares a of the J2 estimation and 48 % for (3. Less than 
2 years of VEX data have a bigger influence than a large interval (more than 30 years) of 
accurate Mars observations. This is especially relevant for the PPN parameter (3 with an 
improvement of about 48% of the accuracy between a determination including only Mars 
data and another one with both Mars and VEX data. In the same time, the improvement 
induced by the addition of Viking data is about 20 % for the J2 and 10% for the PPN 
parameter f3. These figures show the crucial role of the VEX data before the use of future 
data from the ongoing generation of Mercury orbiters. 

2.3. Incremental method and sensitivity estimation 

An original strategy to study the sensitivity of the planetary ephemerides to 32 and 
PPN f3 is to estimate how does differ from INPOP08 an ephemerides built using different 
values for 32 and PPN (3 and fitted on the same set of observations as INPOP08. Such 
differences give an indication on how observations are sensitive to these parameters and 
with what accuracy can we estimate a parameter such as (3. 

To test such sensitivity, we focus our attention on the postfit residuals of the most 
accurate dataset used in INPOP08 adjustement: the Mercury direct range because of 
its sensitivity to general relativity and to the sun 32, VEX, MEX and MGS/MO data 
because of their high accuracy and simulated S/N presented on table [T] and the Jupiter 
Galileo data and Saturn Cassini normal point. These 2 latest data sets are selected 
because they induce a global improvement of the planetary ephemerides in its all and 
especially of the Earth orbit. 
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Table 3. f3 interval in which the residuals stay below the 5% limit (upper table) and 
comparisons between several other estimations found in the literature (bottom table) . In 
the case of (Williams et al. 2009), the value of (3 given here is estimated for 7 = 1. 
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To estimate the sensitivity of these 7 most accurate sets of data used in INPOP08 
adjustment to the variations of values of J2 and PPN /3, we have estimated and plotted 
the ratio S/N defined as: 

fo,o 

where aij is the 1-sigma dispersion of the postfit residuals of an ephemerides based on 
INPOP08 but with values of J2 and PPN /? different from the ones used in INPOP08 
(which are P = 1.0 and J2 = 1.82 x 10"'^) and fitted on aU the INP0PG8 data sets and 
(To.o is the 1-sigma dispersion of the postfit INPOP08 residuals. We have used 9 values 
of J 2 varying from 1.45 x lO"'^ to 3.05 x lO"'^ with a 0.2 step and 24 values of PPN 
(3, building then 192 different ephemerides. The 24 values of f3 are distributed over 2 
windows: a global one based on 12 values of /? varying from 0.997 to 1.003 with a 0.0005 
step (window 1) and from 0.9996 to 1.0004 with a step of 0.0001 (window 2). Results 
presented as the S/N pourcentage, are plotted on figures [2] and ??. 

As one can see on figures the impact of the PPN f3 is not symmetric to /3 = 1. On 
figure [21 one notices also the direct correlation between the S/N obtained with MGS/MO 
and MEX data and the one obtained for VEX. 

One may see on the figure [2] that the S/N of the Jupiter and Saturn data sets 
are sensitive to changes in J2 and PPN /?. The sensitivity of these datasets are not 
crucial for the analysis but they reflect the impact of the use of such observations in 
the improvement of the Earth orbit and then the sensitivity of the Earth orbit to the 
gravity testing. On tabled we have gathered minimum and maximum values of PPN /3 
defining the sensitivity interval of the different datasets. The sensitivity interval is the 
interval of PPN (3 for which the S/N remains below 5%. Values of PPN (3 greater than 
the maximum value given in table [3] or smaller than the minimum value cannot be seen 
as realistic in comparison to modern observations. By considering the figures [2] and ?? 
and the table [3] it appears that the MGS /MO and MEX data provide the most narrow 
interval of sensitivity with 0.99995 < (3 < 1.0002. This interval is in agreement with the 
latest determinations done by (Williams et al., 2009), (Fienga et al 2008) and (Pitjeva 
2006). 
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Fig. 2. Residuals obtained by comparisons between observations and ephemerides esti- 
mated with different values of PPN /? (values given on x-axis of each subframes) and 
different values of Sun J2. 
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3. Secular advances of planetary perihelia 

We are interested here by evaluating if the observations used to fit 1NPOP08 would 
be sensitve to supplementary precession of the planet orbits. Such precessions would 
have be seen as supplementary secular advance of the orbit perihelion unexplained by 
general relativity. To estimate the sensitivity of the modern tracking data, we first fix 
J2=1.8 X 10^^, /? = 1 and 7=1. By fixing the value of the sun 32, we then isolated the 
impact of the secular advance of the perihelion for one given value of 32. 

For each different values of tt, initial conditions of planets are fitted to the INPOP08 
observations and we compare the postfit residuals to the INPOP08 ones. We focused our 
study for the same sets of observations as for the 32,(3 study. As one can on figure [3l the 
behaviour of the obtained S/N (as defined in section [2. 3p is symmetrical to a minimum 
value, this minimal value being centered around n = or not. This symmetry explains 
why in table[3]we gives interval of tt for which the minimum of S/N is obtained. One can 
then compared these values to those published by (Pitjeva 2009). For all planets, except 
Saturn, the values of tt minimizing the residuals are not significantly different from zero. 
One can note that the best constraint on the Earth orbit is given by the Jupiter VLBI 
datasets which gives the narrowest interval of tt. For Saturn, an offset in the minimum of 
the S/N is obtained for the Cassini tracking data sets (—10 ± 8) and the VEX datasets 
(200 ± 160). These estimations lead to two statistically significant and different deter- 
minations of a supplementary precession of the Saturn orbit. By comparisons, (Pitjeva 
2009) value is very close to the one we obtain in considering only the S/N induced on the 
Cassini observations. This result shows how important is the description of the method 
used to evaluate such quantities. 

To test the stability of the estimations and as it is well-known that the asteroids induce 
a global precession of the inner planets perihelia, we operate the same computations with 
small changes in the mass of the ring (20%) and in the Sun 32 (5%) values. The obtained 
variations of the S/N are plotted on figure [3] where the red curves are the results obtained 
with the change in the mass of the ring and the blue curves are the ones deduced from 
the 32. Some changes are noticeable for Viking and Jupiter, however, for Cassini and 
VEX, the minimum are stable. 

The investigation about a statistically significant advance in the Saturn perihelion 
has to be continued in using more Cassini and VEX data. Indeed, a prolongation of 
the interval of time covered by these two datasets will improve the accuracy of the 
estimations. 

4. Does the Pioneer anomaly impact the ephemerides ? 

Since 2002 and the confirmation by several teams of the detection of acceleration anoma- 
lies in the tracking of several spacecrafts, three classes of possible explanations were 
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Table 4. tt intervals minimizing postfit residuals. 
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Fig. 3. Residuals obtained by comparisons between observations and ephemerides esti- 
mated with different values of n 



given; first, the detected acceleration is not really an acceleration but is an avatar or a 
manifestation of a mis-modeling in the Doppler and ranging signals taped by navigation 
teams. Second, the anomaly is a mis-modcling in the orbit of the probe itself induced by 
a technical problem or misunderstandings of the spacecraft techniques. The third cause 
invocated is a generalization of the second by implying a mis-modeling in the dynamics of 
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Fig. 4. Residuals in right ascention and declination of Neptune and Uranus obtained 
with INPOPOS (solution of reference) and a fitted solution including PA of different 
magnitudes: from 8 to 2 xlO^^*^ m.s^^ . The x-axis are years and y-axis is in arcseconds. 



the probe but also of all objects in the solar system and beyond. Thus, if the equivalence 
principal is followed, the equations of motion of the major planets of our solar system 
have also to be modified in the same manner as the spacecraft dynamical equations are. 

We investigate this question by using the INPOPOS planetary ephemerides as a test 
bed for some hypothesis describing the pioneer anomalies. 

A classic description of the pioneer anomalies (PA) is the apperearance of an constant 
acceleration of about S.75 x 10^^*^771.5^^, Sun-oriented after 20 AU (Anderson et al. 
2002a). We then add this constant acceleration in the equations of motions of Uranus, 
Neptune and Pluto. This constitutes the modified ephemerides as noted in the following. 

We have fitted the modified ephemerides to observations usually used to built 
INPOPOS. Residuals obtained after the fit are plotted in Figure |4l The value of the 
acceleration was changed in a way to obtain a minimum value for which the effect in- 
duced by such acceleration becomes detectable in the residuals. As it appears clearly in 
the residuals of Uranus right ascention, a constant acceleration of S x lO^^^m.s^^ added 
to the classical Einstein-Hoffmann equation of motions can not be missed even after the 
fit of the initial conditions of planets. A systematic effect remains especially after 1930. 
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This effect cannot be absorbed by the fit or by the noise of the oid Uranus observations 
made at that time. By changing the value of the acceleration, one sees that the acceler- 
ation must be at least 4 times smaller than the one commonly adopted to be absorbed 
by the residuals. For Neptune and Pluto, the situation is different. For these planets, the 
effect of a constant acceleration is absorbed by the fit, as one can see on figure |4] with 
the postfit and prefit residuals of Neptune. 

5. Conclusions 

Concerning the determination of the PPN parameter /?, an estimation of the planetary 
ephemerides sensitivity to this parameter is done foUwing two methods. Our results show 
that a global fit is needed in order to decorrelate parameters such as PPN /3, Sun 32 and 
the asteroid pertubations. 

We have tested possible detection of anomalous advance of perihela of planets. More 
investigations are needed for the analysis of the perihelion rate of Saturn and more 
observations of Cassini and VEX data are necessary. 

Finally, the results obtained here for the Pioneer Anomaly induce that no constant 
acceleration larger than 1/4 the PA can affect the planets of our solar system. If it was 
so, it would have been detected sooner. In the frame of the equivalence principle, this 
means that no constant acceleration larger than 1/4 the PA can be realistic. 
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Abstract. In this paper, we present several gravity tests made in using the last 
INPOP08 planetary ephemerides. We first propose two methods to estimate the 
PPN parameter /3 and its correlated value, the Sun J2 and we discuss the correla- 
tion between the Sun J2 and the mass of the asteroid ring. We estimate possible 
advance in the planet perihelia. In the end we show that no constant acceleration 
larger than 1/4 the Pioneer anomaly can affect the planets of our solar system. 

1. Introduction 

Thanks to the high precision achieved with the observations deduced from the tracking of 
spacecrafts, it becomes possible to estimate relativistic parameters, for instance mainly 
7 and /?, of the Parametrized Post-Newtonian formalism of General Relativity (Will, 
1993). Nevertheless if 7 plays a role in the equations of motion, it is worth to note that 
light propagation is only sensible to that parameter. PPN 7 can be then estimated with 
high accuracy by light deflection measurements by VLBI (Shapiro et al. 2004), by time 
delay during an interplanetary roundtrip and by Doppler tracking data of a space mission 
(see for instance the Cassini experiment, Bertotti et al. 2003). This is also why, in the 
following, we put 7 = 1 in order to only test the sensitivity of PPN (3 on the perihelion's 
advance of planets. However a physical quantity relative to the Sun, its oblateness J2, 
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plays also a key role in this phenomena. Indeed, usual expression of the advance of 
perihelion is given by (Will 2006) 

_ 277(27 -P + 2)GM,un 37r J2fiLn 

a(l - e2)c2 ^ a2(l - e^y ^ ' 

where G and c arc the ncwtonian gravitational constant and the speed of light in a vac- 
uum, respectively; J2 a-nd Rgun ^rc the Sun oblatcness and the Sun equatorial radius; 
a and e arc the semi-major axis and the eccentricity of the prcccssing planet. One can 
immediately suspect that it is not possible to do a relevant estimation of PPN (3 with- 
out considering Sim J2. Furthermore, it may be not possible to decorrelate safely these 
two quantities with only one planet. With INPOP08 (Fienga et al. 2009), MEX and 
VEX tracking data have lead to an important improvement of Mars and Venus orbits, 
respectively. It is then suitable to take advantage of this new situation by attempting to 
decorrelate these parameters. 

Indeed, the impact of the Mars and VEX observations is not only limited to the 
improvement of the planet dynamics. They also play a role in the determination of pa- 
rameters such as the asteroid masses, the oblatcness of the Sun and the PPN parameter 
/3. The ratio between the uncertainties of the observations and the sensitivity of the 
observed orbit to the GR modifications has been evaluated by dividing the cumulative 
advance of the perihelion over a period of time corresponding to the time span of obser- 
vations by the angle uncertainty of INPOP and presented in table ??. If the amplitude 
of the advance of the perihelion on Venus and Mars orbits is considered for a set of 
observations of equivalent accuracy, Venus data will be seven times more efficient to test 
general relativity and to estimate the sun .12 than Mars. If VEX mission is prolongated 
from 2 years to 4 years and if VLBI observations are done from the tracking of the 
spacecraft with an accuracy of about 1 mas, VEX data will be then as important for the 
PPN testing and Sun J2 estimations as the direct 80G-meter accuracy radar ranging on 
Mercury. Besides, the Mars data are still very important because of the long time span 
of observations of very good quality obtained since the Viking mission in 1978. 

Thanks to the informations brought by the combination of very accurate tracking 
data of spacecraft orbiting different planets, the planetary ephemcrides become then an 
interesting tool for gravity testing. In the following, we give some examples of such tests. 



2. Determination of PPN /3 and sun oblateness J2 

2.1. Correlation between sun J2 and asteroid modeling 

The advance of the perihelion induced by general relativity and sun J2 has an impact 
very similar to the advance induced by the main-belt asteroids on inner planet orbit. In 
INPOP08, a ring was fixed to average the perturbations induced by the main-belt aster- 
oids which cannot have their signal fitted individually on tracking observations. This ring 



